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FOREWORD 


This  report  is  submitted  as  part  of  Contract  AF  33(038)  211;06, 
which  was  initiated  by  the  Office  of  Air  Research.  Work  under  this 
contract  was  be  pun  at  the  Graduate  School  of  Aeronautical  Engineering 
at  Cornell  University  in  March  1951*  The  report  is  one  of  a  series  to 
be  published  as  the  result  of  the  work  carried  out  under  this  contract. 

The  project  is  being  directed  by  Professors  Nicholas  Rott  and 
W.  R.  Sears,  and  the  constant  assistance  of  Professor  Sears  on  this  in¬ 
vestigation  is  especially  acknowledged.  The  investigation  was  suggested 
to  the  author  by  Professor  Arthur  Kantrowitz,  whose  interest  in  Mach 
reflection  originated  in  his  researches  on  "Propagation  of  Expansion  and 
Compression  Waves,"  which  are  partially  supported  by  the  U.S.  Office  of 
Naval  Research. 

Mr.  L.  S.  Wasserman  of  the  Flight  Research  Laboratory,  Wright  Air 
Development  Center,  was  the  project  engineer  of  this  report, authorized 
by  Research  and  Development  Order  No.  465-5-6,  Practical  problems  in 
Aerodynamics. 
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ABSTRACT 


The  configuration  of  non-stationary  Mach  reflection  has  a 
conical  similarity  in  time-space,  since  there  is  no  fundamental  time 
or  space  interval  involved.  This  property  is  easily  shown  (Part  lb), 
and  can  be  used  to  simplify  the  determination  of  various  disturbance 
cjiantities. 

For  this  type  of  reflection  the  strength  of  the  reflected  shock 

depends  on  that  of  the  incident  shock  and  also  on  the  deflection  angle. 

To  the  first  approximation,  this  reflected  wave  is  a  sonic  front*  Thus, 

the  resulting  boundary-value  problem  in  linearized  theory  is  relatively 

1,2,3 

simple,  and  has  been  attacked  by  several  investigators  •  The  ex¬ 
tension  to  second-order  theory  is  discussed  (part  Id)  but  is  not  carried 
out  in  detail. 

Based  on  the  results  obtained  from  linear  theory,  it  is  shown 
how  the  shock  strength  can  actually  be  determined  to  the  second  order 
(Part  II).  The  strength  of  this  reflected  shock  is  found  to  be  of 
second  order,  and  it  vanishes  at  the  triple  point. 

PUBLICATION  REVIEW 

The  publication  of  this  report  does  not  constitute  approval 
by  the  Air  Force  of  the  findings  or  the  conclusions  contained  therein. 

It  is  published  only  for  the  exchange  and  stimulation  of  ideas. 

FOR  THE  COMMANDING  GENERAL; 

B.  vnLLI^ST^olonel,  tJSAF 
-^^^0hief,Flight  Research  Laboratory 
X  Research  Division 
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I.  Mach  Reflection 


(a)  General  Consideration 


Let  us  consider  a  plane  shock  X  propagating  into  still  air 
with  a  constant  velocity  in  the  positive  X  direction;  the  air 
behind  the  shock  has  a  uniform  velocity  TTt  .  Let  us  further  denote 
the  flow  region  ahead  of  the  shock  by  o  ,  and  that  behind  by  /  • 

With  respect  to  a  frame  attached  to  shock  Z  ,  this  is  a 
stationary  flow  field*  The  strength  of  shock  I  is  measured  by  the 

pressure  ratio  =  f**  / p  0  ,  or  by  the  Mach  number  in  flow  / 

ZTS  -T7t 

referred  to  frame  /  ,  i.e.,  t7~  =.  - 

^ '  XT,  T7S 


Since 

_jl  _  2.  IT'CT-  Cd^i) 

/  "  IT's-  / 


V, 


i  ft-  a~xJ 

*  p - —  4,  , 


-  cr1"- 


gVr 

*  <r- 


(i- 


Us  - 


o~  ' 


r° 

/°o 


two  immediate  conclusions  can  be  drawn. 

i)  Since  ^  >/,  i.e.,  0~  <  I  ,  with  respect  to  flow  in  /  ,  shock 
L  always  propagates  with  subsonic  speed,  so  that  any  disturbance  in  / 
catches  up  with  the  shock.  Thus  when  the  shock  hits  a  wall  corner,  the 
appearance  of  Mach  reflection  is  natural  unless  the  turning  angle  $  be 
so  large  that  the  so-called  "regular  reflection"  occurs. 
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ii)  The  Mach  number  of  flow  /  is  Riven  by 

T7,  <r\> 

=  If,  ~  fJV,)  <r 

Thus  the  flow  behind  shock  is  supersonic  < M,>t)  -when 

<r  <  J  (*&>'■->  - 


2.  V 


■  94 


or  ->  z  S' (r^-  .9 

and  subsonic  (  <t  ")  when  or  >  e 6  ,  or  JT  <  ^ , 

■When  such  a  propagating  plane  shock  hits  a  corner,  the  con¬ 
figuration  of  resulting  Mach  reflection  can  be  stated  as  follows: 

There  is  a  triple  point  ~r  ;  part  P 7~of  shock  Z  remains  straight, 
being  unaffected  by  the  signal  from  corner  A  j  part  is  curved,  due 


to  effect  of  deflection  by  A  1  TP  is  the  reflected  wave  signaling 
the  presence  of  A  (in  the  first  approximation,  the  entire  front  is 
sonic  in  subsonic  case,  while  only  part  Tp’  is  sonic  in  supersonic 
case);  and  a  slip  stream  T5  appears,  across  which  entropy  and 
density  (but  not  pressure)  are  discontinuous.  (Fig.  1,  2).  The  whole 
region  propagates  to  the  right  with  velocity  ZJ,  and  at  the 
same  time  expands  radially  about  origin  o  with  speed  of  sound  . 
Inside  region  the  flow  is  non-steady.  However,  by  a  proper 

transformation  of  coordinates  it  can  be  shown  that  there  exists  a 
conical  property,  with  o  as  apex,  and  time  ~t  as  axis  of  the  cone, 
so  that  the  problem  can  be  reduced  to  a  two  dimensional  steady  one. 


The  flow  will  then  be  signified  by  two  parameters  in  its  conical  plane: 
strength  of  the  incident  shock,  and  angle  of  wall  corner. 
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p 


p 


(3)  subsonic  (ii)  supersonic 

An  exact  solution  to  the  problem  is  mathematically  not 
practical,  so  a  linear  solution  is  first  sought.  This  has  been  done 
by  different  authors  under  different  assumptions. 

If  $  is  very  small,  reflected  wave  is  weak.  Vorticity  J~L 
above  the  slipstream  can  be  neglected  (indeed,  SI,  <  Of*?*)),  However, 
unless  incident  shock  is  also  weak,  vorticity  in  region  0  $  is  by 
no  means  negligible. 

(1) 

By  assuming  a  weak  shock  and  small  corner  angle,  Bargmann 

investigated  the  disturbance-velocity  distribution  in  the  whole 

region  'Co  by  introducing  a  velocity  potential.  This,  however, 

necessitates  an  assumption  that  the  shock  be  so  weak  that  the 

existence  of  the  slip  stream  and  the  region  of  rotational  flow  can 

(2) 

be  neglected.  Later,  not  assuming  weak  shock,  Lighthill  showed 

that  the  disturbance  pressure  P  in  whole  region  satisfies  one 

single  conical  equation  in  spite  of  the  presence  of  the  slip  stream 

and  the  rotational  flow,  and  so  he  calculated  the  pressure  distribution 

(3) 

along  both  the  wall  and  Mach  shock.  Recently,  Ludloff  ,  using  the 
Lorentz  transformation,  determined  both  pressure  and  density  field 
in  •  His  analysis  showed  the  slip  stream  lies  along  the  radial 

line  o~T  ,  as  expected,  flow  in  region  o  being  rotational,  while 
that  in  region  o'1"  P  irrotational. 
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Utilizing  the  results  of  previous  authors,  it  should  he  possible 
to  determine  the  disturbance  velocity  a  i\f  in  region  o  f  P  ,  taking 
into  consideration  the  slip  stream,  sc  that  the  incident  shock  need 
not  be  weak,  although  the  deflection  angle  would  still  be  assumed 
small.  This  solution  should  prove  an  immediate  improvement  over 
Bargmam^s. 

Although  we  allow  the  shock  strength  to  be  not  small,  we  shall 
in  the  present  report  still  limit  it  so  that  the  flow  behind,  the  shock 
is  subsonic j  the  case  where  the  after-flow  is  supersonic  will  be  left 
to  further  investigation. 

Essentially,  the  procedure  followed  will  consist  in  taking  the 
uniform  flow  behind  the  straight  shock  as  our  basic  flow,  and  de¬ 
flection  angle  A. }  as  a  perturbation.  Thus  all  the  perturbing  quantities 
will  be  expressed  in  orders  of  4^  •  Based  on  the  results  of  the 
present  investigation,  further  approximations  can  be  carried,  out  using 
a  technique  suggested  by  Lirhthill. 

(b)  Conical  Similarity 

Designate  the  uniform  flow  field  behind  the  propagating 
straight  shock  by  /  ,  the  disturbed  flow  behind  the  curved  Ikch  shock 
and  inside  the  reflected  front  $  by  2  •  Let  velocity,  pressure, 

density  and  entropy  be  denoted  by  9  ,  />.,  /?,  .  (  J  =  /,  2  ), 

i  ' 1  t  c 

and  the  corresponding  non-dimensional  perturbed  quantities  by  ( u,  y),  £ , 

so  that 
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%x.  -  firvi  ,  'J') 

Pi  -  Pi  +  **'?•&/* 

f>2  =  <°,  ( >+  p> 

A*  i  -  1  >  /  Q, 

Then  referred  to  the  fixed  coordinate  system  -with  corner  as  origin. 


the  equation  of  continuity,  equations  of  motion,  and  equation  of 


energy  (isentropy  along  flow  line)  take  on  the  following  forms 


*  ?,  o+*>%2;  +  % 


a-t 

arU 


’Z 


s/°'  0  c?  s'314  '22' \ 

^Jy  *  n  7,  (Jx  t  ~>y)  =■  o 


*  ?,  «-  7 ,  "Zy  +  K' 


~~Z3  * 


air 

y?  ~  "t 

On  introducing  the  transformation 


—  © 


x  - 


3-7'* 


s 


Y_ 

<9.  t 


(1) 


(2) 


we  are  referring  the  flow  to  a  new  reference  frame  which  has  its  origin 
at  o  and  moves  together  with  flow,  and  we  are  contracting  the  whole 
flow  field  by  a  scale  Q.’t  .  This  means  that,  at  any  instant,  the 

t 

sonic  front  that  originated  from  O  in  the  physical  plane 


reduces  to  the  unit  circle  in  the  transformed  plane: 


/  - 


=  / 


(U) 


If  by  this  transformation  we  can  eliminate  the  time  dependence,  then 
our  problem  simplifies  to  a  two-dimensional  steady  one.  That  is,  the 
non-steady  motion  will  have  a  conical  property  in  space-time ;  the  apex 
of  the  cone  being  at  o  ,  the  axes  of  the  cone  being  the  time  ordinate, 
normal  to  the  plane  of  the  motion. 
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Now  since 


3  (  2)  ^  — —  — . 

«?z  =  9T  " 

Jt  ~  ~  (  1?  *  **,  ~  * 

we  have,  on  introducing  the  transformation,  the  following  system  of 
eqiations  of  continuity,  motion  and  energy 

JO  f  -  /sy/>  C  -rTx  *  ) 

_L  a  P 

£>  «  —  /*-  f°  x 

, 

u  -  rz?  "V 


(5) 


-IK  l/. 


where 


JO  f>  -  c/+  M.f> >  (  jz  +  Sj  ) 

h  +  "  A/*  ^  “A  ^ 


=  y 


The  form  of  the  system  of  equations  clearly  shows  that  the  motion  is 
indeed  conical  in  time-space. 

Using  vector  notation  s.  (  k  ,  v-)  t  j X  ~  l**  these  take 


the  form 


(r  -  m. 3  ).v %  =  ^ 

(X  -*>.$)  ■  vP  =  “<■ ^  ? 

<-*-  -  =  “  -+  v  ? 


(6) 


(  V  refers  to  JT  plane  hereafter) 

(c)  Linearization 

Let  us  suppose  that  all  the  perturbation  quantities  and  their 

derivatives  are  small,  i.e. 

k  ,  'S~,  /t  /*  -  &Jc.  <  <  / 
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Then  we  have 


±  ■  V J°  —  M,  V-  ^ 

r  •  rf-  =  */> 

X  '  V  p  =  £ 


(7) 


or 


BP 

'Vfl  _ 

A*/ 

d-x 

+  7  ~ 

2-JL 

J  X 

*  pi}  = 

■V  x 

-%  */ 

3? 

+  #  -0/ 

-z>x 

+  7*4  = 

ar* 

B  X 

<5  - 


^1/ 


(7*) 


In  conform ty  with  the  approximation,  it  may  be  noted  that  the 
reflected  front  Tt  coincides  with  the  sonic  front,  i.e,,  r-  /  ,  and 

slip  stream  5  T  coincides  with  radial  line  6  —  Ow  '  C  ~  ^ • )  • 
By  eliminating  the  undesired  variables,  it  is  easily  obtained, 
i)  Pressure  perturbation 

V-  r.  v  3  =  ( O  r  ?  -  ( r-X7i-t)  (jr-Vf>> 

^  V  > 


i.e. 


=  ( r  'V  r  * ~T)i”  ) 


(&) 


Thus  even  when  vorticity  is  present,  satisfies  a  differential 
equation  which  is  well  known  from  usual  steady  conical  flow  problem.  The 
problem  is  immediately  solved  by  introducing  the  Chaplygin  transformation, 
with  boundary  values  determined  by  consideratioh  of  Rankine-Kugoniot 
conditions.  Accordingly  we  have 


/•  ' 


-  o 


(9) 
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where 


$  « 


/  - 


T 

is  the  Chaplygin  transformation,  ■‘■his  is  the  approach  by  which 

(2) 

Light hill  calculated  the  pressure  distribution. 

(ii)  Density  perturbation,  since 


V  f  -  M.  P 


(10) 


(11) 


so  in  a  region  where  there  is  no  vorticity,  it  is  immediately  obtained 
that  J°  —  /*•  However,  it  should  be  noted  here  that  C*  is  now 

divided  into  two  separate  regions  by  the  slip  stream,  and  this  nice 
simple  relation  holds  only  in  region  °~r  f?*  For  region  °  p  has 

different  boundary  condition  than  f>  and  the  solution  will  accordingly 
be  different.  This  fact  is  clearly  born  out  by  Lu dl off's  calculi tion. 
(iii)  Velocity  perturbation 

or  nr-  =  v-t  , 

Cr-  =  o 

By  elimination,  it  is  easily  obtained  that 

\ 7 l  C c  O  ( ?•  v)  -  ^  o  (12) 

Now  { jr.  &  +  X  )  S  -  ° 

gives  j-  _  cp(g) 

So  cn  introducing  Ghaplygin  transformation  (10),  we  obtain  in  general 

?  =  f(s)  <pce>  (13) 

i.e.,  satisfies  a  Poisson  differential  equation. 
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However,  if  the  flow  is  irrotational,  elimination  again  results 
in  equation  (9). 

r'/  =  l  r.  17*0  f  r.  (n3,) 

or 

Thus,  if  we  confine  our  attention  to  the  vorticity  -free 
region,  ov  p  #  it  is  seen  that  the  problem  simplifies  to  solving  the 
equation  (13’)  subject  tofollowing  boundary  conditions: 

i)  along  wall  o  V  ,  u  -  °  ,  ^  >  ii)  over  T’T  >  i.e., 

r  ■=.  I  ,  u  =.  v-  —  O  •  and  iii)  on  o~r  ,  normal  velocity 

component  vanishes  since  this  is  a  slipstream, 

(d)  Higher  Approximations 


From  above,  it  is  clear  that  with  boundary  conditions  properly 
determined,  all  first  order  quantities  u  ,  \a  ,  p,  can  be  calculated. 
Now  suppose  these  have  been  found,  A  second  approximation  can  be  readily 
carried  out  in  following  manner: 


Putting 


$  -  T 


J  +  3  -f 


where  S  denotes  any  of  m  ,  \s- ,  p ,  p  ,  and  collecting  the  second 


order  terms  in  eqs. 

(6),  we  obtain 

_f  *  f1^-  M>  \ 

A*,  tP  + 

?  njJ 

r  • 

r  •  ^ 

(15) 
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Thus  the  equation  for  the  second-order  term  in  becomes 

*  -  '  -9  'J 


ry";.  (r  *.)  (r  §  * 

-  „  -  ■  "w-v  v'" 

-  M, 


A  fV 

B  t>€^ 

-  [  £  V  v  t^r  +  r)  f'°J  f  *"  ^  J 


-  ( 


r,>  v  <U 


/ 


J> 


(16) 


or 


The  corresponding  boundary  conditions  over  the  wall  can  be  ob¬ 
tained  by  standard  procedure,  those  along  Mach  and  reflected  shocks 
should  be  obtained  by  first  determining  the  second  order  values  on  the 
original  boundary  using  Lighthill 1 s  technique.  Thus  the  work  involved 
in  a  second  order  solution  will  in  general  be  prohibitive,  and  accordingly 
not  attempted.  Only  the  shock  strength  of  reflected  wave,  which  is  of 
physical  interest,  is  investigated  in  next  part  of  this  report. 
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II  Strength  of  Reflected  Shock 


Since  for  small  wall  deflection  angle  &  the  flow  behind 
straight  shock  and  slip  stream  is  irrotational,*  we  are  justified 
to  introduce  a  velocity  potential  $  ,  which,  on  account  of  the 
conical  property  shown  before,  further  reduces  to  the  form 

(  r  *  9,4)  =L  4  T  *  f  r  r'  g)  (!) 

Vp,  _  X— 

where  we  have  put  <3,  =.  ~jcT  >  r  ~  <?f  +  • 

Now  denoting  partial  derivatives  by  subscripts,  the  radial  and 
tangential  velocity  components  u  ,  \/~  can  be  written 


£  r*  -  -fy  <3 , 

v-  — 

and  since 

§*$&  ~  T  £ 9 

c£t)r*  « 

f  £-e)r  +  rt 

** 

t 

% 

il 

JX  J*  a  1 

So  we  have 
i) 


1  r  («'+  ^  ^  *  rZ 

y-  f  fir--) //-  rA-  •*  i  (2) 

ii)  «%  f  -  +  S  Sxf 

f  v*{  =  <3,y  ^  i-i)  ah- 


*  n  <  o(*)y)  Bargmann,  p.  16 
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as  j 

i  -  4  *'f 

v-t- 

a,  r-r 

®  -  f  +rr 

K'’ 

a,r  f  _ 

-  T 

r  {kG  ) 

i.e. 

jrfstrel 

a‘  2/ 

— v  ^  f  - 

r'frr  +  *-  C~ 

r-fr/>r  + 

4,%  1 

■*  f  7 C  Vv,r  * 

Atf 

rJ 

+  r*~  & 

rv 

Eliminating  from  i)  and  ii) 

_  rxfr,  -  ir/r  -frr  +  fg  -  h.f<yfrO 

*  f,Vrr  *  7? +  ■k-Sr/e/'rO 


$  ,  -  r  l-<  r-->)  C -f  *  {  r-J*]  <-  r-f*  -  IS-'j} 

+( -fer^ni  A/  •frOs  -f&&  l  ~  ° 


Referring  to  coordinates  fixed  with  respect  to  shock,  the  Rankine- 

Hugoniot  condition  can  be  written 

.  „  U) 

i )  7  —  continuous 


ii) 


_i  f  —  -  *«>! 

r-*i  L  * 
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•where  ^  1  and  ^  <t,}  denote  velocity  components  tangential  and 
normal  to  the  shock  respectively,  and  )  denotes  the  increase  of  a 
quantity  across  the  shock  front. 

Now  since  <£>  exists,  and  flow  is  steacty-  in  the  Crs  6 ^  plane, 
let  the  reflected  shock  front  be  r*  =  (&)  ,  then 

«  ^  a,  ( -  1  +  -fr)  *  -  y  fe 
^  (t)  /&  -h  v1  (9)  z_  T  f  _ 


M  J~  t-t-  r  *v  ']'x' 

(•>)  _  ~  Cfr-V>  +  ’ 

/  '  +  r '*■*?' 


■x  7'  fe 


i)  gives 


[  *  i  're>  - 


i.e. , 


ii)  -  -d/i-  ~r-xr'*fr  .  _  I - — .  ^ 

A.  *?  - -  - -■-  ■  - -  *»,  —  -J  K»  Jr 


0  J7*  r-'f''- 

*  ,  ax  _  »-  r  ■  ->K»lJ-rJr>  ,„,j 

ZT 1  '  b  ~  r*<  1  %<"/«?-  1 


Therefore  we  have 


tuic  no  iictvo 

^  v  7  -  fr  *•  *  *  *7  jB  7 

^  =  -  /T,  IzFTJTTT,  ^ 

=  .  J-  [  ,  .  -  f.+  (i-o-f'l  J 

f  <£«*/)  -  ( 

=  _i_  [ar-.U-f-  rfr)  -  */r+  *-<?-'>] 


Since  /,  -  //*  -  ° 
in  linear  theory. 


—  /  and  •/*,  q  ~  A  (O')  9t 
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r-h 


—  -f  r/« 


a  (8  )u/7 ^  n.  -*  f 


Following  light hill  by  putting 

r  *  R+  *~f9)  =  *  +  +  ^(8)  +  ■  ■ 

-f  ^  /,  I  /?,  d)  +  +••• 

in  (If)  and  condition  on  &  -  (  ,  that  the  linear  coefficients  in 

(li )  must  vanish;  and  noting  that  although  since 


^  -*  ^  “  r'/9)*n 


terms  and  /gg  may  produce  terms  in  ,  they  are  however 

all  not  of  dimension  1;  we  have 


-2  Rr,  - (£- * 


on  f2  -  / 


i.e,  r,  (Qy  s  O 


/>»*/)  *//*  //e/2  =  /’/'*  O  /?///?  ^ 

'fxn  *  fi.  OQ  -  "T~/  M*  A  '~(8) 

l-xr^  -  rr-'M-St-  /**)  a 

'“x  ^  /*  V<?J 
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The  shock  condition  at  ^  are:  Since  o 

ahead  of  shock,  so  directly  behind  it 

r  r  _  *  *n  -o 

(12) 

Shock  cones  in  range  fi<  t ,  y  <  ✓  -t  r\.  ■+•  <j  ■+ 
so  7  ~ 1  °i  dimension  2,  7  *■!'*'  7+  +  f  i  j  °  ^  7 '  r*  ■ 

In  condition  /*=- o  on  r=-  ,  terms  /» ,  •  .can  be 

expanded  by  Taylor  expansion  (at  R  =  l)  using  fist*  f*]x,  -r*.}  ■*  *  • 

Terms  of  dimension  1  and  2  give 

=  -r^A-fO)  (13) 

r-x  =.  v*; 

To  express  the  second  condition  of  (2  ),  /*  is  needed  when 
fist-*  -  rLJ  +  f  ^ j  -  O J> 

The  largest  terms  here  are  of  dimension  2,  and  these  come  from 


,/,-a?  and  from  (-fxa) 


n.  x.  i 


so  they  are 


a*,  ^  A  (8)  J  r-^-  ■* 


by  (12)  and  (13) 


m,  xs4  ( e ) 


c  r^-**)  v  /  _ 

Putting  7x  a  “  /w*  f\  fGJ  ,  th 


this  becomes 


/>/£>  /TvffO  -  A70J>  =  -*■  *7<9>>  -  /4W 
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The  solution  of  (15)  gives 


K  (&)  ~ 


A  (&}  >  o 
A  C&>  <  o 


(16) 


For  the  case  of  compression,  A  ( 7>  a  ,  the  shock  is  thus 
located  at 


r  -  t  -+ 


/  6 


(i?) 


New  the  pressure  change  across  the  shock  is 

fif*  =  pjS<rf'-f>  * 

Since  the  shock  strength  is  given  by 


n 


- V  ~  ^  *V/V/>  (X0) 


/»,  ClAn) 


where 


A  (B> 


-  ,J tzJLL—  \ 

r->  ,(  ) 


~  3  n 


on  linear  theory 
in  Chaplygin  plane. 


and  as  given  by  Ref.  2  (59)  is 


*t  ^  /  i?'/  lL  _ j- 

it  at  triple  point  (  -x  a  /?,  /  -  00  is  o  ,so 

strength  of  reflected  shock  vanishes  at  the  triple  point. 
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CONCLUSIONS 


The  conical  property  in  time-space  of  diffraction  of 
propagating  plane  shock  is  discussed.  Second  order  solution  for 
pressure  j?  has  been  suggested,  but  since  this  requires  more  informa¬ 
tion  regarding  n  ,  and  the  second-order  boundary  conditions  than 
is  presently  available,  actual  carrying  out  of  the  solution  is  not 
attempted. 

The  strength  of  reflected  shock,  however,  has  been  investigated; 
and  it  is  found  to  be  generally  of  order  /$  .In  particular,  this 
vanishes  at  the  triple  point.  The  immediate  conclusion  is  that  while 
the  strength  is  of  second  order  in  general,  it  can  be  at  most  of  third 
order  at  triple  point  in  particular. 
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APPEND H 


The  Reflected-Shock  Angle  at  the  Triple  Point 


A  slight  extension  of  this  investigation  discloses 
that  at  the  triple  point  where  the  incident,  reflected  and 
Mach  shocks  meet,  not  only  is  the  reflected  shock  of  vanishing 
strength  hat  also  is  it  locally  tangent  to  the  acoustic  circle# 
The  latter  statement  can  be  proved  as  follows: 

Denoting  by  (//  the  angle  included  between  the 
reflected  shock  and  acoustic  circle,  then 


dr 
Kot  0 


where  r  is  given  by  Eq.  (17)# 

Now  the  mapping  from  the  S  ~  P  Q  ‘ 


plane 

2 


to  the  2,  -  X,  +  d  % ,  plane,  in  Light  hill's  paper  ,  is 
everywhere  analytic,  so  /  /  and  its  derivatives  are 

finite#  Even  though  the  mapping  from  the  physical  conical 
plane  f  -e.  ‘  to  the  5  plane  is  singular  on  the  unit 

circle,  the  following  relations  hold: 


•O  * 


pao 


''Op  y /  -IT  1 


So  it  is  immediately  seen  that  this  singular  behavior  will  not 

affect  the  determination  of  */”  «  furthermore,  since 

a  _  _s.  _  ~z>  /  d?/  / 

/°de  ~at  ~  off! 


we  eventually  get 


'-d  r  _  /,  oL/ c/2ii  “f_A 1 


and  the  conclusion  follows  at  once# 
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